Evidence for Negative Stiffness of QCD Strings 
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QCD strings are color-electric flux tubes between quarks with a finite thickness determined by 
the dimensionally transmuted coupling constant, and thus with a finite curvature stiffness. Con- 
trary to an earlier rigid-string model by Polyakov and Kleinert, and motivated by the properties of 
a magnetic flux tubes in type-II superconductors, we put forward the hypothesis that QCD strings 
, have a negative stiffness. We set up a new string model with this property and show that it is free 

^\ • of the three principal problems of rigid strings — particle states with negative norm, nonexistence 

0^ I of a lowest-energy state, and wrong high-temperature behavior of string tension — thus making it 

a better candidate for a string description of quark forces than previous models. 
■ PACS number(s); 11.17.-|-y, 12.38.Aw, 12.4Q.Aa 

<: 
00 

m 



o 
o 



X 



Color-electric flux tubes between quarks have a thickness of the order of the dimensionally transmuted coupling 
constant. They should therefore display a finite resistance to a change in extrinsic curvature. To account for this 
property, a higher-derivative stiffness term was added to the Nambu-Goto (NG) action of an infinitely thin string 
The resulting Polyakov-Kleinert (PK) string seemed to render a more relevant description of gluon forces in 
K*" QCD than the former NG string. Particularly attractive was the fact that a PK model containing only a stiffness 
C3 term is asymptotically free at short distances and can generate the string tension spontaneously For other 
properties see [Q. 

Unfortunately, however, the PK model has also serious consistency problems. First of all, there exists an unphysical 
ghost pole in the propagator |^ . The pole is generated by the second derivatives with respect to the string coordinates 
\^ in the action. Second, the PK model shares with all higher-derivative theories the lack of a lowest-energy state 
^\ , Third, the squared string tension of the PK model has an unphysical imaginary part at high temperatures or low 
• /3 = l/T, even though this model unlike the NG string displays the correct power behavior for /3 ^ consistent 
I with asymptotic freedom (except, however, for an imaginary factor i) [^. 

It is the purpose of this letter to propose a new string model in which the problems of the PK string are absent 
' while the attractive features are preserved. This new model has a negative extrinsic curvature stiffness. It was 
inspired by properties of magnetic flux tubes in type-II superconductors and analogous properties of Nielsen-Olesen 
vortices in relativistic gauge models. In the London limit, these have a nontrivial energy spectrum, from which several 
calculations have deduced a negative stiffness [p^|-p^. By analogy, we hypothesize that the stiffness QCD strings is 
also negative. The action we propose for a description of such a string is 



A = I d'i^g'^ Ax^ z il5vJ^ D,xx , (1) 

where x^(^) with A ~ 0, 1 are the string coordinates in a d-dimensional euclidean spacetime parametrized 

by * = Oi 1j and 9ij — diX^djXx is the induced metric on the world surface of the string with being its inverse 
and g = det{gij). Covariant differentiation with respect to ^* is denoted by Di, and = DiD^ is Laplace-Beltrami 
operator. The physics of this model is governed by two mass scales M and ^, and a dimensionless constant c. 

Due to nonzero radius of thickness the interaction between surface elements of the string has a nonlocal character, 
just like in a magnetic flux tube. In momentum space, the quadratic part of the action (j^) gives rise to the free 
propagator 

G{k') = . (2) 



For small fc, this has an expansion 



G{e)='-±^ + ..., (3) 
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with the mass parameter = (c — The stiffness a is the inverse coefhcient of in the action and has the 

negative value a = —h^jM"^. 

The propagator contains only a single pole at fc^ = 0. This is to be contrasted with the PK model where 

corresponds to the positive stiffness a = Is? jM"^ and contains an unphysical pole with negative residue at fc^ = — A^. 

The inverse propagator (H) is plotted in Fig. ^ for sample parameters /i and c. The fluctuations with large momenta 
are more violent than in the PK string with the opposite curvature stiffness, but less violent than in a NG string with 
reduced tension (c — l)M^/c, which is eventually approached. 

The above string is inspired by the properties of a vortex line in a type-II superconductor [^3| , which has a diameter 
of the order of the penetration depth governing the parameter 1/A of the model. The inverse mass corresponds 
to the coherence length which goes to zero in the London limit. 

In QCD, the parameters of the model have the following physical origin. As long as there are no quarks, QCD 
gives rise to a string with only one length scale which enters into both the string tension and the thickness of the 
flux tube. Thus it determines the two parameters M and A with some ratio between them, which is a consequence 
of the detailed dynamics of the gluon system. As quarks are included, this ratio is modified and an additional length 
scale arises at which the vacuum polarization of the gluon propagator becomes important and changes the stiffness 
properties of the flux tube. This length scale is proportional to some average over the inverse quark masses, and 
accounted for by the model parameter 

For practical calculations, it is convenient to approximate by (|^), and cut all momentum integrals off at k'^ = jj? . 
To have a physical penetration depth, we must take A to be smaller than the cutoff /x, or 1 < c < 2. 

The action associated with the approximate Green function (|^) is 

^1 - Ul^ I d?i^g'^ D,x^ ^ _ ^^^^^ D,x^ . (5) 




FIG. 1. The inverse propagators of the model (|l|) (solid curve) and of the approximation (g) (dashed lines) for fi = 10 and c = 1.2. 
The dotted curve shows the behavior in a PK string with the opposite value of the curvature stiffness, the dash-dotted curve corresponds 
to a NG string with reduced tension (c — l)Af^/c. We see that the large-fe fluctuations are much more violent than in the PK string, but 
less violent than in a NG string, which is eventually approached. 

We shall now derive the high-temperature limit of the total string tension, which can be done exactly in the limit 
d oo. As usual in this limit, we make the metric field gij an independent fluctuating field, and force it to be equal 
to the induced metric diX^djXx by means of Lagrange multiphers A*-' , adding to (^) a term 

A2 = ^M^ I d^^^y' {d.x^djxx - 9,j) . (6) 

It is straight-forward to calculate the partition function of the string at a finite temperature T anchored to two static 
quarks separated by a large distance R. The world sheet of the string is periodic in the imaginary-time direction with 
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a period f3 
sheet as 



T ^. Choosing coordinates with running from to R, and from to /?, we parametrize the world 



with x°'{^), a = 2, 3, . . . d— 1, describing the transverse fluctuations which occur quadraticafly in the extended action 
Ai+ A2- These can be integrated out producing a trace log with a factor d—2 which in the limit d ^ oo suppresses 
the fluctuations of gij and A*-' , so that the saddle point approximation provides us with an exact partition function. 
At the saddle point, we may assume diagonal forms for the metric and the Lagrange multipliers: 



with the constants pi and A^. The resulting effective action ^cff =-^cff + A^°ff^ consists of a mean-field term 



(7) 
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(8) 



which also appears in NG and PK string models, and a one- loop contribution 
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(9) 



Here the temporal components ^o ^re summed over all thermal Matsubara frequencies go = (9o)n = 2?™/ (3yfpQ, n = 
0,±1,... . 

The effective string tension is defined by M"^^ — F//3, where F — Acs/R is the free energy per unit length evaluated 
at the extremal values of pi and A^. At low temperatures, /3 3> 1/A and the one- loop term (H) is very small compared 
to the mean-field term so that the extremum of AcS hes at po = pi = 1, Ao = Ai = 0, where F — j31vp . At 
high temperatures, on the other hand, /? <C 1/A and the one- loop term becomes essential. In this limit, it can 
easily be evaluated analytically. 

First, however, let us make an instructive observation: If we assume for a moment that the cutoff p? be artificially 
small, so that + q\ ^ A^, we can neglect the denominator of the first term in Eq. (||) and see that the reduced 
one-loop action (||) takes the NG form. By the same calculation as in Rcf. |l4j we then find at high temperatures the 
squared free energy per unit length 



F' 



{d-2)n 



(10) 



This is in bad qualitative disagreement with the behavior derived from QCD in the limit of large- TV by Polchinski 



^qcd(/9) 



(11) 



In contrast to this, the present model has the correct high-temperature behavior. Since the cutoff /i^ is, of course, 
much larger than A^, the negative curvature stiffness a = — A^/M^ can take effect. In the high-temperature limit 
f3 <C 1/A, the first term in the logarithm (^ is small compared to the second-derivative term, and one- loop action 
(^ reduces to 



/l'°°P ~ - - R 




(12) 



The first term comes only from zero Matsubara frequency, in the large fjl^ -limit, the second contains the full spectral 
sum in (^), but approximated by the leading second-derivative term. Adding to (|l^) the mean-field term (|^), the 
resulting high-temperature effective action AcS has a saddle point at 
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Po = 2 - Ai , Pi = 2 + Ao . ^ . , (13) 

Ao + Ai Ao + Ai 

where Aq and Ai are determined from 

^|^y^(Ao + AO = l, ^^iPo'-^)=^4^Po'- (14) 
In terms of Ai, the associated squared free energy per unit length reads simply 

= M^P'^iX, + if . (15) 
Pi 



The deconfinement temperature is determined by the vanishing of 1 / pi , which yields 

-1/2 

(16) 



Mtot(O) {d-2) 



Svr 



Here Mtot(O) = M \/v + 1 is the total string tension of at infinite 6eia, and the parameter v is related to the stiffness 
a as follows 



(d-2)|a| 



The parameter 6(|a|) is an unique real positive root of the polynomial equation of 9th order which exists for a 
reasonable range of the mass parameter A > 8.25il/. Numerical studies show that the deconfinement temperature 
( p^ ) depends only weakly on the stiffness parameter a and lies in our model always around Td/Mtot(0) ~ 0.932. This 
is somewhat larger than the estimates derived from rigid strings 

Equations ( p^ ) and ( p^ determine Ai via a polynomial of 14th order. In the above range of the mass parameter, 
A > 8.25Af , there are real roots for all /3. At high-temperatures, these have the asymptotic behavior 

A, = const. X |const.|>2; ^ ""^^^ ■ d^) 

Inserted into (p^), they yield the diagonal elements of the metric. The resulting high-temperature behavior of ( p5| ) is, 
therefore, 

'K\d-2f (Ai - 2)(Ai + 1)^ 

36/32 ^^(x^ - 1)2 ■ ^ ' 



This this agrees in sign and /3~2_bghavior with the QCD result (11), if we neglect the weak /3-dependence of the 
running coupling constant g{ji) in (pT|). 

Note that the two strings are not immediately comparable since the result (^l|) was obtained in the N oo -limit 
where the QCD string becomes infinitely thin and has therefore no curvature stiffness. At large but finite N, however, 
it does acquires a finite thickness, as we know from Monte Carlo simulations, and the behavior ( pi] ) can be modified 
at most by small terms of order 1/A^. In this regime the agreement between our result (|l|) and the QCD result (|l^) 
is definitely significant. 

We consider this agreement as an important evidence for a negative stiffness of hadronic strings, which thus behave 
very similar to vortex lines in type-II superconductors. Our conclusion should be verified in lattice gauge simulations. 
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